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The status of simulations using the non-perturbatively O(a) improved Wilson action in the quenched approxi- 
mation is reviewed. The impact of non-perturbative improvement on the hadronic mass spectrum and the size of 
residual lattice artefacts in spectral quantities and decay constants are assessed. 



1. INTRODUCTION 

One of the most important systematic effects 
in lattice simulations of QCD is the finiteness of 
the lattice spacing a. It is well known that physi- 
cal observables computed using the Wilson action 
are subject to corrections of order a, which can be 
rather large. In order to obtain reliable results in 
the continuum limit it is desirable to reduce lat- 
tice artefacts, either through the Symanzik im- 
provement programme , or by employing a 
renormalisation group approach |3[^] . Recently 
the ALPHA Collaboration has carried out the 
Symanzik improvement programme to leading or- 
der through a non-perturbative determination of 
the 0(a) improved fermion action and isospin 
currents |^-||] . This approach should lead to the 
complete removal of lattice artefacts of order a 
in spectral quantities and matrix elements of lo- 
cal currents, so that the remaining cutoff effects 
are of order a 2 . The non-perturbatively im- 
proved action has already been employed in a 
number of simulations in the quenched approxi- 
mation Here, we assess the impact of non- 
perturbative improvement on the calculation of 
the mass spectrum and decay constants in the 
light hadron sector, analyse the scaling behaviour 
and estimate the size of residual lattice artefacts. 

The general expression of the 0(a) improved 
fermion action reads 

+ Csw T iivF,j, v {x)il){x), (1) 

where Sp is the (unimproved) Wilson fermion ac- 
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tion and c sw is an improvement coefficient. The 
improved and renormalised axial and vector cur- 
rents are defined as [jl3| 

(A a )« = Z A (1 + b A am q ){Al + c A ad^P a }, (2) 
(F R )» = Z v (l + 6 v am q ){u; + Cyad^}, (3) 

where Z A ,Zy are the renormalisation factors of 
the respective currents, and 6a, ca, by and cy 
are further improvement coefficients. The unim- 
proved currents and densities are defined as 

P a = ? 7 5^V>, TZ v =Sfa lu ,$1>, (4) 

where r a are the Pauli matrices acting in flavour 
space. The normalisations Z A ,Zy and improve- 
ment coefficients c sw ,ca,cv,&v have been deter- 
mined non-perturbatively for bare couplings go 
in the range < go < 1 1|-§| ■ Results and pro- 
posals for non-perturbative determinations of b A 
at go — 1 have been reported |l4]-|il| . Further- 
more, all of the above improvement coefficients 
have been calculated in perturbation theory to 
one- loop order 0,^8| ■ 

The parameters of the simulations discussed 
here are listed in Table In addition to the non- 
perturbative value of c sw , the QCDSF collabora- 
tion have also used unimproved Wilson fermions 
for a direct comparison §. UKQCD have also 
used the tadpole improved value of c sw on the 
same set of configurations and at (3 = 5.7 |p9| , po| , 
where a non-perturbatively determined value of 
c sw is not available. Table shows that all col- 
laborations have used the same two values of (3. 
This implies that one cannot test as yet whether 
non-perturbative improvement indeed leads to a 
scaling behaviour of physical observables which 
is consistent with 0(a 2 ) corrections. However, if 
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Table 1 



Simulations using non-perturbative improve- 
ment. The number of exceptional configurations 
discarded from the ensemble is shown in brackets. 
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Figure 1. Edinburgh plot. 



one assumes this to be the leading scaling be- 
haviour after improvement, one can still estimate 
the size of residual lattice artefacts at a given 
value of a. 

2. SPECTRAL QUANTITIES 

In Fig.[j] we show the Edinburgh plot. In ad- 
dition to the non-perturbatively improved data, 
we also display the tadpole improved results by 
UKQCD at (3 = 5.7 and the data by GF11 using 
the unimproved action at (3 — 5.7, 5.93, 6.17 pl|. 
One observes that improvement yields consis- 
tently lower values for mjv/m p compared to the 
unimproved action. The most dramatic effect is 
observed at /3 = 5.7 when one compares the unim- 
proved results (full triangles) to the tadpole im- 
proved ones (full circles). In fact, it seems that 
the mass behaviour obtained using tadpole im- 
provement at (3 — 5.7 is indistinguishable from 
the non-perturbatively improved action. How- 
ever, as we shall see later, the residual lattice 
artefacts in the tadpole improved data at (3 — 5.7 
are still large. Thus, the Edinburgh plot disguises 
rather than exposes lattice artefacts and should 
therefore not be used to draw conclusions about 
the scaling behaviour. 

It is a well-known fact that lattice simulations 
using unimproved actions fail to reproduce the 
experimentally observed behaviour of the vector- 
pseudoscalar mass splitting, i.e. my — fripg ~ 
const, which holds up to the mass of the charm 
quark. This is usually ascribed to the influence 
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Figure 2. UKQCD data for the vector-pseudo- 
scalar mass splitting in units of m 2 K . Asterisks 
denote the experimentally observed values. 



of lattice artefacts in the computation of these 
splittings. Data for m v — m| s by UKQCD |To[ ] 
and QCDSF Q using the non-perturbatively im- 
proved action show that the splittings are close to 
the experimental values for the (p, tt) and (K* , K) 
systems (see Fig.|| for a plot of the UKQCD 
data). Furthermore, by comparing the results at 
[3 — 6.0 and 6.2 one observes that the dependence 
on the lattice spacing is small. However, for quark 
masses above m s trange the two collaborations see 
a slight downward trend in the data, so that one 
can expect that the D* — D splitting is still not 
reproduced correctly. 

We now discuss the chiral limit and the critical 
value of the hopping paramter, k c . Usually n c is 
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defined at the point where the pseudoscalar mass 
vanishes, rnps = 0. In accordance with the quark 
mass behaviour of mp s implied by the PCAC re- 
lation, one can determine k c from a linear fit to 

(amps) 2 = a 2 B\ (± - = a 2 B m q . (5) 

Both QCDSF and UKQCD have reported that 
this linear ansatz in l/n results in poor fits with 
large correlated x 2 /dof and have therefore re- 
sorted to using model functions which also con- 
tain quadratic terms in However, in order to 
be consistent with 0(a) improvement, the quark 
mass m q in eq. (|^) should be replaced by 

fhq = m q (l + fe m am q ), (6) 

where b m = —J — 0.0962<?q in one-loop perturba- 
tion theory ]18|. This modification of the fitting 
ansatz in the determination of k c has so far not 
been used. 

Another method defines k c at the point where 
the quark mass defined through the PCAC rela- 
tion vanishes, i.e. mpcAC = 0. Here it is impor- 
tant to realise that even after non-perturbative 
improvement chiral symmetry is only approxi- 
mately restored, so that 

<9 M {A^x) + c A ad il P{x)} 

= 2m? CAC P{x) + 0{a 2 ). (7) 

Therefore, the values of k c determined by requir- 
ing either mp s = or mpcAC = will differ by 
terms of order a 2 . The compilation of results in 
Table]^ shows that in the range of (3 under study 
the difference in n c using either method is of the 
order of 10~ 4 and thus statistically significant. 

We now analyse the scaling behaviour of the 
mass of the vector meson by comparing the ap- 
proach to the continuum limit for unimproved 
and improved actions. To this end, we note that 
the continuum limit should be studied for con- 
stant physical volume, so that finite-volume ef- 
fects do not distort the scaling behaviour. Since 
chiral extrapolations are poorly understood, we 
will use the available lattice data for the vector 
and pseudoscalar masses omv,amps as well as 
the string tension a^fa in order to interpolate 
the dimensionless ratio my/^/a to mps/v^ = 



Table 2 



Values for k c determined at mp s = and at 
"IpcAC = 0. 
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0.13589(2) 0.135861(19) 
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ALPHA 
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6.0 


0.135196(14) 




6.2 


0.135795(13) 


0.135802(6) 



1.125. For sfd = 440 MeV this implies m PS = 
rriK = 495 MeV. If mp s and my are linear func- 
tions of the quark masses then one expects that 
m v ^ m K * = 892 MeV or my/V^ - 2.027. In 
Fig.|| we have plotted the interpolated data for 
my/y^a versus a\fa for unimproved ]2l|] and im- 
proved actions (ll] Jl^,^] ■ For the unimproved ac- 
tion one observes that my ~ rriK* is satisfied, 
but only after the extrapolation to the contin- 
uum limit. In contrast, the data obtained us- 
ing the non-perturbatively or tadpole improved 
actions at (3 = 6.0,6.2 show very little depen- 
dence on a and are rather close to my ~ nix*- 
This, however, can no longer be claimed for the 
tadpole improved data at /3 = 5.7, where in- 
stead one observes large residual lattice artefacts. 
From the slope of the linear fit to my / yfa for 
the unimproved data, we infer the size of residual 
lattice artefacts at a ~ 0.1 fm to be (12 ± 1)%. 
For the non-perturbatively improved data we es- 
timate the leading corrections of 0(a 2 ) to be only 
~ 2%. 

We can now turn the tables and ask how well 
various prescriptions to fix c sw satisfy the "phe- 
nomenological" improvement condition 

my/y/a — const (8) 

Besides the data obtained using tadpole and non- 
perturbative estimates for c sw , one can use fur- 
ther spectrum data obtained using c sw = |2l| 
and c sw = 1 p2fl . In Fig.|] we plot my/^ff at 
mps = rriK versus c sw . From the plot one infers 
that the condition (||) is satisfied within statisti- 
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Figure 3. Scaling behaviour of the vector mass 
at tops = rriK- The dashed line indicates the 
extrapolated result using the GF11 data. 



cal errors if 

> tadpole _ 



> 6.0, 



(9) 



where u denotes the average link variable. Thus, 
with the present statistical accuracy, already the 
tadpole improved estimate of c sw leads to a large 
reduction of O(a) effects, provided the lattice 
spacing is not too large. Fig.|] also shows that 
in order to satisfy ^ at — 5.7 one would have 
to go to much larger values of c sw than those im- 
plied by the tadpole prescription. This, however, 
appears impossible due to the appearance of ex- 
ceptional configurations Q . 

3. DECAY CONSTANTS 

From the definitions of the renormalised and 
improved currents, eq. (|l]), one sees that in or- 
der to obtain O(a) improved matrix elements of 
axial and vector currents one requires knowledge 
of the improvement coefficients ca, 6a, cy, by and 
the renormalisation factors Za and Zy. In this 
section we will focus on the pseudoscalar decay 
constant (/ps)r, which we write as 



(/ps)r = Z A (l + b A am q ){f^ + c A g^}, (10) 

where the unrenormalised matrix elements are 
parametrised as 



mp S /W = <0|A 4 (0)|PS) 
g£i = <0|P(0)|PS). 



(11) 
(12) 




Figure 4. my/^/a as a function of c sw . The sets 
of points from right to left have been obtained 
using non-perturbative, tadpole, tree-level and no 
improvement. Data at — 5.7 are shown only for 
unimproved and tadpole improved actions. The 
dashed line shows the extrapolated GF11 data. 



A major advantage of non-perturbative improve- 
ment is that the renormalisation factor Z A has 
been determined with an accuracy of around 1%. 
At = 6.0 and 6.2 one finds § 



Zk = 



0.791(9), 
0.807(8), 



= 6.0 
= 6.2 



(13) 



Note that the error on Z \ must be combined with 



f(0) 



(0) 

CAflps- 



the statistical error on / pg 

The improvement coefficient ca has also been 
determined non-perturbatively H : 



ca 



— 0.083, 
— 0.037, 



= 6.0 
= 6.2 



(14) 



The effect of ca on the decay constant can be 



studied by comparing /pg with fp§ +CAffpg • Here 
the contribution from CAffpg leads to a decrease 
in the pseudoscalar decay constant of <~ 4% at 
= 6.2 and even ~ 15% at = 6.0. This effect 
is particularly pronounced for < 6.2. 

The improvement coefficient 6a is relevant for 
the determination of the kaon decay constant fx . 
However, unlike by, the coefficient 6a has so far 
not been determined non-perturbatively. In one- 
loop perturbation theory one finds f\M 



(0) 



,(0) 



6 A = 1 + 0.1522^ + 0^). 



(15) 
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In order to study the influence of 6a, one can 
evaluate (/ps)r around the strange quark mass 
for different choices of 6a- Here we compare 

• 6 A = 

• b A = 1 + 0.1522gg 

• b A = by 

The seemingly ad hoc choice of 6 A = by (here 
we use the non-perturbative determination of by) 
is motivated by the observation that the one- 
loop coefficients for 6a and 6y are approximately 
equal |Q . By applying different choices of 6a to 
the analysis of the UKQCD data around mx, one 
finds that 6a leads to an increase in (/ps)r of at 
most ~ 2% at /3 = 6.2 and - 3% at /3 = 6.0, 
which is fairly small. Since the choice 6a = by 
gives essentially the same mass behaviour com- 
pared to choosing 6a = 1 + O.I522<7q, one con- 
cludes that perturbative estimates of 6a are quite 
acceptable for quark masses up to and around 

^strange • 

We now analyse the scaling behaviour of 
(/ps)r- In Fig. | we plot (/ P s)Rr versus 
(mps?~o) 2 , where tq is the hadronic radius defined 
in 1 23 1 . Data for r /a were taken fro m p4| . If lat- 
tice effects are small the data in Fig. |5j should lie 
on a universal curve. The results for Jk^q com- 
puted at P = 6.0 and 6.2 show a slight dependence 
on the lattice spacing. In order to study residual 
lattice artefacts, we employ a similar procedure 
as in the case of the vector mass and extrapolate 
(/ps)rH) to (rriKro) 2 - Fig.|| shows the resulting 
values of fx^o as a function of (a/ro) 2 . There is 
good agreement between the data from all three 
collaborations, and in principle their results could 
be combined. Furthermore, it appears that a lin- 
ear extrapolation in a 2 yields a continuum result 
which is compatible with the experimentally ob- 
served value (although there is a priori no reason 
why the quenched approximation should repro- 
duce the measured result). In contrast, the au- 
thors of |p5| have found that, for the unimproved 
action, the continuum value of Jk is significantly 
lower than the experimental result. 

From the slope in (a/ro) 2 one estimates that 
residual lattice artefacts in /k^q amount to ca. 




Figure 5. The pseudoscalar decay constant as a 
function of TOpg in units of ro for j3 = 6.0 (open 
symbols) and (3 = 6.2 (full symbols). Asterisks 
denote the experimental values for and fx us- 
ing ro = 0.5 fm. 
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Figure 6. JkTq plotted versus (a/ro) . 
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10% at a ~ 0.1 fm. Compared to the previously 
discussed case of the vector meson, this is a fairly 
large correction. Given the substantial contribu- 
tion of CAffpg to (/ps)r at (3 = 6.0, one can ask 
whether ca has a large effect on the scaling be- 
haviour. If one formulates a similar improvement 
condition for ca as in eq. (||), for instance 



fxro = const. 



(16) 



one can study how well it is satisfied for different 
choices of ca- Using the UKQCD data it turns 
out that the above condition "favours" smaller 



6 



values of ca at the lower end of the /3-range. How- 
ever, at this stage one should not jump to con- 
clusions before a more thorough investigation of 
improvement conditions has been performed. 

4. CONCLUSIONS 

First results from simulations using the non- 
perturbatively improved Wilson action and cur- 
rents in the quenched approximation show that 
lattice artefacts are drastically reduced. Unlike 
the case of unimproved actions, the results for 
the vector meson mass are practically indepen- 
dent of a at (3 = 6.0 and 6.2, so that residual a 2 
effects are around 2% at a ~ 0.1 fm. A real test of 
the scaling behaviour is still lacking and will only 
become possible when more and different values 
of the lattice spacing are studied. 

The analysis of data for the pseudoscalar de- 
cay constant has shown that a non-perturbative 
determination of 6a is required for quark masses 
above ra s trange- The improvement coefficient ca 
has a large influence on the scaling behaviour, 
which motivates further investigation. On the 
whole, it appears that non-perturbative improve- 
ment leads to better agreement between the con- 
tinuum result for fx and the experimental value. 

The systematic nature of non-perturbative im- 
provement makes it easily applicable to other sit- 
uations: results for c sw computed for two flavours 
of dynamical quarks have been reported |2C| . 
Furthermore, the formalism has been applied 
to quenched QCD with an improved gauge ac- 
tion §7|. 

1 thank P. Rowland for help with the UKQCD 
data and also M. Guagnelli, T. Mendes and G. 
Schierholz for communicating their data prior to 
publication. I wish to thank M. Liischer and R. 
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